Self-Energy Correction to the Two-Photon Decay Width in Hydrogenlike Atoms 
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We investigate the gauge invariance of the leading logarithmic radiative correction to the two-photon decay 
width in hydrogenlike atoms. It is shown that an effective treatment of the correction using a Lamb-shift "poten- 
tial" leads to equivalent results in both the length as well as the velocity gauges provided all relevant correction 
terms are taken into account. Specifically, the relevant radiative corrections are related to the energies that enter 
into the propagator denominators, to the Hamiltonian, to the wave functions, and to the energy conservation 
condition that holds between the two photons; the form of all of these effects is different in the two gauges, 
but the final result is shown to be gauge invariant, as it should be. Although the actual calculation only in- 
volves integrations over nonrelativistic hydrogenic Green functions, the derivation of the leading logarithmic 
correction can be regarded as slightly more complex than that of other typical logarithmic terms. The domi- 
nant radiative correction to the 2S two-photon decay width is found to be —2.020 536 ( a /w) (Za) 2 ln[(Za)~ 2 ] 
in units of the leading nonrelativistic expression. This result is in agreement with a length-gauge calculation 
[S. G. Karshenboim and V. G. Ivanov, e-print physics/9702027 |, where the coefficient was given as —2.025(1). 

PACS numbers: 12.20.Ds, 31.30.Jv, 06.20.Jr, 31.15.-p 



O 
Of 



> ■ 

in : 
o . 
o ■ 

-xt- : 
o . 

o : 
• i— i 

Or 



X 



I. INTRODUCTION 

The two-photon decay of the metastable 25 level in 
atomic hydrogen and hydrogenlike systems is a rather in- 
triguing physical phenomenon; it was first investigated by M. 
Goppert-Mayer a long time ago (jj 01 • The transition involv- 
ing two quanta limits the lifetime of the metastable 25 reso- 
nance, at least for low and medium nuclear charge numbers Z. 
By contrast, the highly suppressed magnetic dipole transition 
to the 15 ground state has a negligible influence on the decay 
width [3]. In this article, we evaluate the dominant self-energy 
radiative correction to the two-photonprocess. We recall here 
the known leading-order result m^HaLl 



r = 8.229 Z b s~ L = 1.310 Z b Hz . 



(1) 



For ionized helium (Z = 2), rather accurate experimental ver- 
ifications of this result exist I8,l9l ll0ll . Due to its metastability, 
the 2S level in hydrogenlike systems is one of the most accu- 
rately defined resonances found in nature. Indeed, it is this 
very property — the small natural linewidth — which has made 
possible the high-resolution two-photon spectroscopy of the 
15-25 transition O El El El El El- 

The fully relativistic quantum electrodynamic formalism is 
intricate when applied to bound-state problems 11171 M 8l ll 9t 
12011 . but it is often possible to gain a rather good understand- 
ing of QED radiative corrections to a particular process if one 
uses a simplified, NRQED Lagrangian that contains effective 
operators which then lead to the perturbations that have to be 
evaluated (see e.g. ll2lll22il '). Of course, the main difficulty of 
any bound-state calculation, which is the separation of the two 
energy scales (scale of binding energy and the energy/mass 
scale of the free particles), persists in the effective approach. 
It is necessary to also specify cutoff prescriptions; the arti- 
ficially introduced scale-se para tion parameters then cancel at 
the end of the calculation 12 ll l22l l23l 12411 . Elucidating dis- 
cussion of the latter point can be found in |25] Ch. 123] and 
in EfJCh. 11.4 on p. 493]. 



Within nonrelativistic quantum electrodynamics (also re- 
ferred to as NRQED, see 1191 [20]), one has the choice be- 
tween two different forms of the interaction Hamiltonian: 
the "length" (Yennie) and the "velocity" (Coulomb) gauges. 
There are certain intriguing issues involved with the gauge 
invariance in the dynamical nonrelativistic atom-light interac- 
tion. Indeed, in order to prove gauge invariance for dynami- 
cal processes, it is in many cases necessary to carefully con- 
sider the gauge transformation of the atomic wave function in 
addition to the transformation of the fields. Otherwise, non- 
gauge invariant results are obtained off resonance I27ll28ll29ll . 
In the current situation of radiative corrections to the two- 
photon decay width, we will show that it is possible to ig- 
nore the transformation of the wave function: the two-photon 
decay width, including the radiative corrections, is invariant 
under a "hybrid" gauge transformation |28] which involves 
only the fields, but ignores the gauge transformation of the 
wave function. In general, the choice of the gauge and the in- 
terpretation of physical operators have to be considered very 
carefully in time-dependent problems (see 113 (A p. 268] and 
Refs. iBElEi). 

The gauge invariance of the two-photon decay rate and 
of the radiative corrections to this effect can be regarded as 
slightly problematic, partly because the integration over the 
photon energy is restricted to a finite interval. By contrast, the 
gauge invariance of the low-energy part of the one-loop self- 
energy shift, in an effective NRQED treatment, holds only be- 
cause one may drop terms whose divergence, for large photon 
frequency, is stronger than logarithmic |23, see Eq. (3.4) ff.]; 
in this case gauge invariance would be violated over finite in- 
tervals of the virtual photon frequency. It has been one of the 
main motivations for the current paper to study related ques- 
tions. 

This article is organized as follows: In Sec. [II] the lead- 
ing nonrelativistic contribution to the two-photon decay rate 
is discussed, together with its relation to the NRQED two- 
photon self-energy. In Sec.|III] the leading logarithmic radia- 
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tive correction to the two-photon decay rate is formulated, the 
discussion is based on a perturbation with an effective poten- 
tial. Explicit expressions are derived in Sees. IIVI and Ivl for 
the length and velocity gauges, respectively. Gauge invari- 
ance is proven in Sec. lVIl Numerical results are presented in 
Sec. IVHI Conclusions are drawn in Sec. IVIIII All derivations 
are presented in some detail, for the sake of transparency. 



II. LEADING-ORDER TWO-PHOTON DECAY RATE 

The decay width of a bound system may be understood nat- 
urally as the imaginary part of the self energy 1 3 1 ] . Indeed, 



the (negative) imaginary part of the self-energy is just T/2, 
where T is the decay width. We discuss the derivation of the 
two-photon width based on this concept, within nonrelativistic 
quantum electrodynamics |32]. 



The formulation of the two-loop self-energy problem 
within the context of nonrelativistic quantum electrodynam- 
ics (NRQED) has been discussed in 13 311 . We denote by p> 
the Cartesian components of the momentum operator p = 
— iV. The expression for the two-loop self-energy shift 
reads | HH 



A-Enrqed — — 

1 / ; 



/ 2a \ 2 f'l f'2 ( I , 1 . 1 , 1 . 

^ / du)\ Wi / dcj 2 <*>2 { ( p p p p 

\3Trm 2 ) Jo Jo l\ H - E + u>i H - E + wi + w 2 H - E + uj 2 



p p p 



2 \ H — E + Ui H-B+U1 + U2 H - E + u 

1 / , 1 j 1 j 1 

p p 



2 \ H - E + uj 2 H - E + u>i + u 2 H - E + uj 2 



H - E + 0)i \H — E ) H - E + lu 2 

(p' ~r, — ^r~, p'*) (p j ( ~Tr — lr~, 1 p'} _ t; { " 



2 \ H-E + uJi I \ \H-E + lu 2 J I 2 \ H - E + uj 2 / \ \ H - E + uj 

_ m I i 1 1 i\ _ ro / ( 1 i\ _ m / i 1 i\ 1 

m \P h - E + ui H - E + u> 2 P / ui+u 2 V H - E + uj 2 P / ui+u 2 V H - E + u>i P / J 



(2) 



All of the matrix elements are evaluated on the reference state 
\(f>), for which the nonrelativistic Schrodinger wave function 
is employed. The expression for the two-photon decay width 
[Eq. @ below] now follows in a natural way as the imagi- 
nary part generated by the sum of the first three terms in curly 
brackets in Eq. (0. Specifically, the poles are generated upon 
u>2 -integration by the propagator 



1 15) reads 



To 



4 a 2 
27 ~ 



H — E + Ui + Ll>2 



T?l 



-f E' - E + ui + uj 2 



(3) 



at LJ2 = E — E' — u>\. Alternatively, this condition may be 
expressed as E — E' = u>i + lo^, and represents the energy 
conservation condition for the two-photon decay. The imagi- 
nary part generated by the first three terms in curly brackets of 
the energy shift (0 is thus seen to yield the two-photon decay 
width IB5II. 

In view of the above discussion, and in agreement with 
Shapiro and Breit [4, Eq. (3)], the nonrelativistic expression 
for the two-photon decay width Tq in the case \<fi) — |25) and 



1 



H -E + UJ2 

j \ 

H — E + lui 



(4) 



where u>2 = w max — u)\ and cj max = E — E' is the maximum 
energy that any of the two photons may have. When compar- 
ing this expression to Eq. (2) of 1 36], it should be noted that 
the quantity y ibid, represents a scaled photon energy. The 
Einstein summation convention is used throughout this arti- 
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cle. Note the following identity [27, 3 



IV. LENGTH GAUGE 



P' 



1 



P' 



m H — E + uj\ Tfi 



+ 



p< 



1 



777 H — E + lo 2 rn 



-U>\ U>2 771 



1 



H - E + lo 1 



H -E 



■ u 2 



(5) 



which is valid at exact resonance lui+uj 2 = E—E', This iden- 
tity permits a reformulation of the problem in the velocity- 
gauge as opposed to the length-gauge form. 



III. RADIATIVE CORRECTIONS 

We consider a hydrogenlike atom and employ natural units 
with h = eo = c = 1. In order to analyze the radiative cor- 
rection to the two-photon decay width, one could write down 
all Feynman diagrams which contribute to the process, and 
start evaluating them. However, a much more economical un- 
derstanding into the problem can be gained by considering an 
approach inspired by effective field theory, or nonrelativistic 
quantum electrodynamics I19l 12011 . in which the leading effect 
due to radiative photons is described by an effective Lamb- 
shift potential Sill 



<5VLamb = ~ a {Za) \n[(Za) ) 



(6) 



In this work we will consider a "standard normalized pertur- 
bative local potential" i40ll 



SV = 



ir(Za) 



S^{r). 



(7) 



(8) 



which is related to <5VLamb by a simple prefactor, 
SV Lanih = \ - \n[(Za)- 2 } SV. 

3 7T 

The corrections to the Hamiltonian, to the energy and to the 
wavefunction, incurred by the perturbative potential Q, read 
as follows, 



E -» 

SE = 

H -> 

10) - 

\H) = 



E + 5E, 

{4>\sv\4>) , 

H + SV, 
10) + 1^} , 
1 



8V\(j>) . 



(9a) 
(9b) 
(9c) 
(9d) 

(9e) 



E-H 

The standard potential Q leads to a "normalized" energy shift 
with unit prefactors, 

{Za) rn 



SE(nS) = 



(10) 



According to (0}, the two-photon decay rate To of the 
metastable 25 state is given by 



A 



where we use the definition 



A- — — 

~ 27 7T ' 

as well as u> 2 = E 2 s - Ei S - u>i and w max 
The quantity ( is given by 



where 



IS 



(2 = (IS 



C = Ci + c 2 



1 



(11) 

(12) 

E 2 s — Eis- 
(13a) 



H - E 2S + uJ\ 
1 

H — Eis — <jJ\ 



25 



25 



(13b) 



(13c) 



The perturbation l|9} leads to the following replacements, 
which include the first-order corrections to the various quan- 
tities that are relevant to the 25 decay width, 



E 2 s 
1 15) 

|25) 

^2 



E 1S + SE 1S , 8E 18 = (1S\5V\IS) , 
E 2S + 6E 2S ., 8E 2S = (2S\5V\2S) ., 



\1S) 



125} 



Eis — El 
1 



SV 



SV 



15 



25 



E 2S -H / 

lu 2 + 8u} 2 , 5w 2 = SE 2 s - SEis ■ (14) 



The latter correction ensures that a perturbed energy conser- 
vation condition is fulfilled, 

oji + lo 2 + 5lj 2 = E 2 s - E 1S + (SE 2 s - SEis) , (15) 

i.e. that the two photon frequencies add up to the perturbed 
transition frequency. 

The first-order self-energy correction <5r to the two-photon 
decay rate may be expressed as 



sr 

B 



^max ^max 

2 J duji ujfuj 2 (SC + 3 Suj 2 J dwi w\ 



3 ,2 z-2 
^2 Sp 



where the correction 8C, is the sum of six terms, 

6 
3 = 1 



(16) 



(17) 
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to be defined as follows, and the second term on the right-hand 
side of d!6t is due to perturbed energy conservation condition. 
The quantity B may be inferred from ©, and l|8} as 

B = --^H{Za)-*]. (18) 



The terms 8Q\ and 6(2 are related to energy perturbations to 
the matrix elements, 



5(i = IS 



1 



, H — E2S + 
5( 2 = (IS\5V\1S) (15 



25 ) (2S\SV\2S) 



1 



H — E\s — uii 

whereas the terms 6(^3 4 5 6 are perturbations to the initial- and final-state wave functions 



25 



(19a) 
(19b) 



6(3 = (IS 
5U = (lS 
S( 5 = (is 
S( 6 = (is 

V. VELOCITY GAUGE 



1 



H — E2S + ^1 \ E2S — EL 



1 



H - Eis —w\ V ^2S _ H 



5V 



5V 



25 



25 



— H J H — E2S + u>i 



1 



1 



1 



25 , 



25 



(19c) 
(19d) 
(19e) 
(19f) 



Gauge invariance of the leading-order decay-rate [see 
Eqs. and £TJ] 



We now discuss the evaluation of radiative corrections in 
the velocity gauge, where the interaction Hamiltonian is given 

by 

rt A A 2 

H[ nt = -e^ + e^. (20) 
to Imr 

According to @ and Q, the leading-order decay rate in the 
velocity gauge is 



r = n 



(23) 



immediately follows from Eq. (|5j; this equation may be 
rewritten in a compact form as 



(24) 



r' 
1 

.4 



dwi uii 0J2 £ , 



(21) 



where A is defined in d!2i . lo 2 = E 2 s — Eig—ui and w n 
£■25 ^ Eis- The quantity £ is the sum of two terms, 



(22a) 



where 



6 
6 



15 
15 



P' 



1 



F 



m H — E2S + u>i m 



F 



1 



F 



m H — Eis — to\ tti 



25 
25 



(22b) 



Equation 1241 may be proven easily by repeated application of 
the commutator relation(s) 

V - = i [H - E 2S + wi, a:*] = i [H - £is - wi, a:*] . (25) 

TO 

Now the first-order correction to the two-photon decay rate, 
in the velocity gauge, is 



5V 

-g- = 2 / dcji uji UJ2 £ 5£ + Slj2 I dwi uj\ £f , (26) 



where the prime denotes the velocity-gauge form of the cor- 
rection and B is defined in ( I18> . We desire to show that 
6T = ST'. 
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The correction <5£ finds a natural representation as the sum In analogy to < !19al and dl9bl i. 6£\ and S& are energy pertur- 
of eight terms, bations, 



(27) 



V 



S£i = (IS 
5& = (IS\6V\ 15) ( 15 



1 



m \ H — E2S + u\ I m 



P' 



25 ) (25 1 S V I 25} , 



2S 



m \ H — E\s — u>i I m 



(28a) 



(28b) 



The terms S^s 4 5,6 are perturbations to the initial- and final-state wave functions, 



= 

^4 = 
<%S = 
% = 



15' 



IS' 



IS 



15 



P' 



1 



P' 



m H — E2S + lui m \ E2S — H 



1 



P' 



1 



V 



1 



m H — Eig — u)\ m \-B2s — -ff 



<5F 



SV 



2S 



2S 



Eis — H I m H — E2S + ^1 m 



SV 



Eis — H I m H — E\s — u>\ m 



25 , 



25 



(28c) 



(28d) 



(28e) 



(28f) 



Finally, S£r,8 are due to the seagull term, 



S& = 15 

m \ 



S& 



1 



-A 15 
m \ 



SV 



E2S — H 
1 



SV 



Exs-H 



25 , 



25 



(28g) 



(28h) 



VI. PROOF OF GAUGE INVARIANCE 



Here, we merely present the results of the analysis carried out in detail in App.|A] Indeed, using Eqs. JAlat — dAlffc . as well 
as jA2at and jA2bL we obtain the compact relation 



<5£ = —uji UJ2 SC — Slo2 oj\ C ■ 



(29) 
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In view of this relation, we can rewrite dl6l and i26\ using ( I29K 



ST' 



2 / du>x Ll>i Ll>2 (— U>1 U>2 C) [ — ^1 W2 <5C — &^2 0J± C] + &^2 / cL>l C 2 




2 / C^C+ (2 + l)6u 2 J duJxulujlC 2 

o o 

ST 
~B ■ 



(30) 



This proves the gauge invariance ST = ST' of the logarithmic radiative corrections to the two-photon decay rate of the metastable 
25 state in hydrogenlike systems. The gauge invariance of the lea ding-order decay rate (Tq = T' Q ) has been indicated in Eg. d23l>. 



VII. NUMERICAL RESULTS 

(Leading order.) We recall that, according to 0, the well- 
known leading-order nonrelativistic effect To is of the order 
of a 2 (Za) 6 . The result for the two-photon decay width of 
the metastable 25 state is 

r = 0.001318 222 a 2 (Za) 6 m. (31) 

This translates into 

T = 8.229 351 997 ZS" 1 (32a) 
= 1.309 742 049 Z 6 Hz. (32b) 

(Radiative correction.) In view of Eqs. (0 and (|8j, the lead- 
ing logarithmic radiative correction ST is of the order of 



ST ~ a 3 (Za) 8 ln[(Za)- 2 }m. 



(33) 



i.e. of relative order a (Za) 2 \n[(Za)~ 2 } with respect to To. 
In the length gauge, the relevant expression for ST can be 
found in Eq. dl 61 . [For clarity, we would like to indicate that 
the correction S0J2 occurring in the expression (I16> is defined 
in (I14> . the quantity £ can be found in Jl 31 . and the terms SQ 
(i = 1, . . . , 6) are defined in Eq. dl9>.l In the velocity gauge, 
the relevant expression for ST' can be found in ( 1261 . with the 
S£i (i = 1) • • • , 8) being defined in Eq. d28i . 

According to dl6l and i26\ . both <5r as well as ST' find 
a natural representation as the sum of two terms, the first of 
which summarizes the perturbations to the matrix elements, 
and the second is a consequence of the perturbed energy 
conservation condition for the transition. Gauge invariance 
ST = ST' has been shown in Sec. lVIl yet it is instructive to ob- 
serve that there are indeed considerable cancellations among 
the two contributions to (5r and ST'. Specifically, we have 
from the first and the second terms on the right-hand sides of 



dl 61 and ( I26> . respectively, 

— = (29.542 - 31.562) - (Za) 2 ln[(Za)- 2 ] , (34) 

r 7r 

— = (8.500 - 10.521) - (Za) 2 \n[(Za)- 2 } . (35) 



[The cancellations appear to be typical for radiative correc- 
tions to decay rates; this has recently been observed in con- 
nection with radiative corrections to the one-photon decay of 
P states l4lll .l 

The final result for the leading logarithmic correction to the 
decay width of the metastable 25 state is 



ST 



<5T' 

— = -2.020 536- (Za) 2 \n[(ZaY 

1 71" 



(36) 



The calculation of S( and 5£ involves expressions analogous 
to those encountered in ll40tl . In ll3^1 . the coefficient has been 
given as —2.025(1), which is in agreement with the current 
calculation. [There is a misprint in the overall sign of the cor- 
rection as given in the abstract of 13611 : one should follow the 
sign indicated in Eq. (8) ibid.] 

The result (I36i . converted to Hertz and/or inverse seconds, 
reads 

.ST = -3.273 x 10~ 7 Z 8 ln[137 2 Z- 2 ] Hz (37a) 



-2.057 x 10" 13 Z s ln[137 2 Z~ 2 } s" 1 



(37b) 



For low Z, the highly suppressed Ml one-photon decay 
25 — > 15 is numerically smaller than the radiative correction 
d37l to the two-photon decay (see Refs. I^I^MI^E^ll ) 
because it lacks the large logarithm: 



Ml 



2.496 x lQ- b Z 



6 7 10 -1 



(38) 



All results indicated in this article for T and ST relate to the 
metastable 25 state; however the approach may easily be gen- 
eralized to the two-photon decay of other states. 
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VIII. CONCLUSIONS 

In the current investigation, the derivation of the leading 
radiative correction to the two-photon decay width of the 
metastable 25 state in hydrogenlike atoms has been based on 
the effective "radiative potential" discussed in Sec. [TO] It 
has been shown that the gauge invariance of the corrections 
holds due to the interplay of corrections to the transition ma- 
trix elements on the one hand and corrections due to perturbed 
energy conservation conditions on the other hand [first and 
second terms on the right-hand sides of ( 1161 and i2bt , respec- 
tively]. The corrections to the transition matrix elements are 
again divided into corrections to the wave function (these were 
referred to as the J-terms in 1 36]), and to the energies that en- 
ter into the propagator denominators, which were termed /in 
the length-gauge calculation I3al . The length- and velocity- 
gauge forms of the correction are discussed in Secs. llVl andlVl 
The gauge invariance of the radiative correction holds (even) 
on the level of the effective treatment as implied by the ra- 
diative potential l|6), as shown in Sec. IVII All derivations are 
presented in some detail, for clarity and transparency. The 
numerical evaluation in Sec. lVIII follows immediately. 

There are two more results of the current paper, probably of 
rather minor importance, which should only briefly be men- 
tioned: first of all, the relativistic result (without radiative cor- 
rections) for the decay rate at Z = 1 has previously been 
indicated as 8.229 s" 1 E3 H |H 0, whereas in jj, the 
(nonrelativistic) result has been indicated as 8.226 ± 0.001 in- 
verse seconds. The current investigation [Eq. J32H confirms 
that the discrepancy has been due to a certain overestimation 
of the numerical accuracy in the early nonrelativistic calcula- 
tion 1 4], not due to a conceivable large relativistic shift. Sec- 
ond, the discussion in Sec. |ll| clarifies that the concept of a 
decay width as an imaginary part of a self-energy 13111 gen- 
eralizes to the two-loop self-energy shift, in which case the 
imaginary part gives rise to the two-photon decay width. 

The leading-order nonrelativistic contribution to the two- 
photon decay width is of the order of a 2 (Za) 6 mc 2 (see 
Sec. QJ. The self-energy radiative correction to the two- 
photon decay is of the order of a 3 (Za) s \n[(Za)~ 2 ] mc 2 , 
as discussed in Sec. IVIII with explicit results indicated in 



Eqs. d36l > and J37l> . It would be interesting to evaluate also 
the constant term of relative order a (Za) 2 . This term supple- 
ments the logarithm evaluated here which is of relative order 
a(Za) 2 ln[(Za)~ 2 ]. According to our experience, in bound- 
state calculations, the nonlogarithmic, constant term has an 
opposite sign as compared to the leading logarithm, and its 
magnitude is two or three times larger than the coefficient of 
the logarithm. This is true for radiative corrections 14111 as 
well as Lamb-shift effects JsHBl. 

One should note a rather general interest in various intrigu- 
ing details related to the two-photon decay process, which are 
not restricted to the search for conceivable parity admixtures 
to the 25 state (see e.g. |53, 54]). Although accurate mea- 
surements of integrated decay rates are difficult 1 55], there is 
some hope that in low-Z and middle-Z ionic systems, exper- 
iments will eventually profit from the possibilities offered by 
electron-beam ion traps, especially when combined with con- 
ceivable x-ray lasers that could be used in order to excite the 
trapped ions into the metastable states. 

Finally, we recall that accurate measurements of the two- 
photon decay width test the 25 state for parity-violating 2P- 
admixtures and can therefore be used as a test for a conceiv- 
able electron or nuclear (electric) dipole moment or for inter- 
actions via "anapole" or "pseudocharge" currents I56ll57ll58il . 
One particularly interesting investigation on hydrogenlike 
Ar 17+ , with an elucidating discussion of the issues related to 
parity admixtures, has been given in 1 59] . 
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APPENDIX A: RELATIONS AMONG MATRIX ELEMENTS 



In this appendix, we present in detail the relations needed for the proof of the identity J29I . For 6£i as defined in j28at . we 
have 



15 



V 



1 



to \ H — E 2 s + wi 

+ (uJ2 - UJl) ( 15 



P_ 

rn 
1 



25 



-U>\ UJ2 ( 15 



1 



H - E\s + 



H - E 2 s + wi 
2S) + (lS\x i x i \2S) . 



25 



(Ala) 



We notice the term <5£i emerge on the right-hand side [see Eq. dl9aH . The corresponding relation for S£ 2 reads 



15 



p' 



1 



P' 



m \ H — E\s — u>i J m 
1 



25 ) = -U1U2 ( 15 



+ {lo 1 - u 2 ) ( 15 



For 6£ 3 , the following relation is useful, 



H — Eis — oji 
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For 5£ 4 , we have 
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The term S£,$ may be reformulated according to 
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Finally, we have for S£ 6 
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However, the relations (IAla> — JAlfl are not yet sufficient in order to proceed with the proof of gauge invariance. We also need 
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We notice the (negative of) the seagull terms (|28g| and J28hl > eme rge. 
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